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ABSTRACT 



Context. We have performed observations of water maser emission towards a sample of low-mass protostars, in order to 
investigate the properties of jets associated with the earliest stages of star formation and their interaction with the surrounding 
medium. 

Aims. The main aim is to measure the absolute positions and proper motions of the H2O spots in order to investigate the 
kinematics of the region from where the jet is launched. 

Methods. We imaged the protostars in the nearby region NGC 1333-IRAS 4 in the water maser line at 22.2 GHz by using the 
VLBA in phase-reference mode at the milliarcsecond scale over four epochs, spaced by one month to measure proper motions. 
Results. Two protostars (A2 and B) were detected in a highly variable H2O maser emission, with an active phase shorter 
than four weeks. The H2O maps allow us to trace the fast jet driven by the B protostar: we observed both the red- and 
blue-shifted lobes very close to the protostar, < 35 AU, moving away with projected velocities of ~ 10-50 km s _1 . The 
comparison with the molecular outflow observed at larger scale suggests a jet precession with a 18' yr _1 rate. By measur- 
ing the positional spread of the H2O spots we estimate a jet width of ~ 2 AU at a distance of ~ 12 AU from the driving protostar. 
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1. Introduction 

The formation of a star is accompanied by a period of 
strong mass ejection in which two lobes of supersonic gas 
move away from the newly born object forming bipolar 
jets, strongly interacting with the natal cloud. The launch- 
ing process of jets is still not well known, nor is the pre- 
cise launch region, i.e. whether the jet originates from the 
interface between the star's magneto sphere and disk or 
from a wide range of disk radii (e.g. IShang et al. 2007t 
Pudritz et al J 120071 ). Most investigations of the jet origin 
so far have been made in atomic microjets from optically 
visible T Tauri stars, thanks t o the high angular resolution 



pered by the high attenuation associated with these high 
density parts of the cloud. 

In this context, water (H2O) masers at 22 GHz repre- 
sent a unique tool to investigate the mechanism of jet for- 
mation and collimation in the earliest star forming phases. 
On the one hand, single-dish surveys show for Class ob- 
jects a definitely higher detection rat e than that found fo r 



the more evolved Class I ones (e.g. iFuruva et alj 120031 ). 



availa ble in the optical range (|Rav et alj|2007t I Bally et al 



120071) . As a next step, it is crucial to investigate the prop- 
erties of jets driven by pr otostars, the so-called Class 
objects ( Andre et al. 1990), which are still hidden in the 
innermost part of the star forming region. However, ob- 
servations of these regions in the optical and even in the 
IR spectral windows are not straightforward, being ham- 



On the other hand, images taken with the VLA interfer- 
ometer have confirmed a close association between wa- 
ter masers and protostars, finding the H2O spots concen- 
trated within severa l hundred AU o f the central star (e.g . 
Terebev et all 119891 ICherninl [l99l iMeehan et all Il998 



Furuva et alj |2005). 



However, there are only a few high angular reso- 
lution studies available, using the Very Long Baseline 
Interferometry (VLBI) , which is needed to trace the inner- 
most part of the jet/disk system at the su b-AU-scales, in 
the cl osest vicinity of th e excit ing object. IClaussen et al 
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(|l998l ) and lFuruva et all (2000) observed respectively the 
HH212 and S106 FIR protostars, tracing with the H 2 
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maser spots the curved bo w shocks though t to tr avel 
along a jet. More recently, iMoscadelli et al. (2006) re- 
ported multi-epoch VLBI observations of water masers 
towards a sample of low-mass protostars in the Serpens 
region a nd towards RNO 15-FI R. For the Serpens Class 
objects. IMoscadelli et al. (2006) found high-velocity spots 
located perpendicular to the jet axis: the authors suggest 
that such maser emission originates at the base of the 
jet, tracing the interaction between the jet itself and the 
accretion disk. For RNO 15-FIR, water masers are dis- 
tributed along a line, suggesting their association with a 
collimated flow, in a greement with the results obtained by 
Marvel et al] (|2008l ) for the NGC 1333-IRAS 4B region. 



Table 1. Coordinates used as phase centers at VLB A 



Thus, the origin of H2O masers in low- mass protostars is 
far from being fully understood, calling for further obser- 
vations at AU-scales. 

Here, we report the results of a multi-epoch H2O maser 
survey towards the NGC 1333-IRAS 4 cluster of proto- 
stars, performed with 0.28 AU resolution and in phase- 
reference mode. Our main aim is twofold: (i) to measure 
the absolute positions and proper motions of the maser 
spots in order to study the kinematics of the disk/jet sys- 
tems, (ii) to provide new observations of low-mass proto- 
stars allowing a step ahead in the discussion of the origin 
of jets. 

2. NGC 1333-IRAS 4 

NGC 1333 is a well-known star forming region containing 
a large number of protostars and located at a distance of 
235 pc , according to the recent observations Hirota et alj 
(j2008l ) performed with VERA (VLBI Exploration of Radio 
Astrometry). The IRAS 4 region is split into three star 
forming sites, called 4A, 4B, and 4C, which can be iden- 
tified in the continuum, from cm- to submm-wavelengths 
(ISandell et alJll99lt iLefloch et all 1 19981: iRodriguez et al ' 



1999). Hi gh angular resolutions cm- and m m-observations 
( Loonev et al. 200Ci Reipurth et al. 2002f ) have revealed 



that IRAS 4A is a binary system composed of Al and A2 
and is driving a highl y collimated outflow in the northeast- 
southwest direction (jChoil 120051 ). IRAS4B is also associ- 
ated wi th a bipolar outflow located in the North-South di- 
rection (|Choill200ll ). IRAS 4C has been rev ealed as a water 



mase r source by single-dish observations (jHaschick et al 



1980| ). althoug h it was not detected by later interfere 



metric survey s ( Rodriguez et al. 20021 : Furuva et al. 2003; 



Park fe Choill2007[) . All these characteristics make IRAS 4 



an excellent laboratory in which to study multiple star 
formation. 

The IRAS 4 star forming reg ion has been observe d in 
the 22 GHz H 2 maser line bv IPark fc Choi (|2007t ) us- 
ing the VLA and a resolution of ~ 80 mas. These au- 
thors found no emission towards the Al, whereas they 
observed six maser spots close to A2. Since the A2 
maser group shows velocit ies close to the systemic ve- 
locity, IPark fc Choi ( 2007 ) interpreted such H 2 spots 
as a tracer of a cir cumstellar disk. On the other hand, 
iMarvel et all (|2008l ) reported VLB A observations of H 2 
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IRAS4A1+A2 


03:29:10.456 


+31:13:32.05 


IRAS4B 


03:29:11.960 


+31:13:08.06 


IRAS4C 


03:29:13.530 


+31:13:58.07 



masers towards NGC 1333-IRAS 4B, finding two groups 
of maser spots expanding away from each other at veloc- 
ities of ~78 kms" 1 along a direction different from that 
of the lar ge scale outflow, alr eady traced with fews H 2 
sp ots bv IPark fc Choi (|2007l ). The disagreement found 
by IMarvel et alj ( 2005 ) calls for further analysis in or- 
der to clarify which are the driving sources and whether 
NGC 133 3-IRAS 4B has a c ompanion as has been pro- 
posed by IMarvel et all (|2008l ). 



3. Observations 

Using the VLB A network, we observed H 2 maser emis- 
sion at 22.2 GHz in phase-reference mode over four epochs 
spaced by about four weeks during 2003 (Apr 1st, May 
14th, June 11th, July 07th). The duration of each obser- 
vation was 10 hours in total. 

As the four protostars (IRAS4A1, 4A2, 4B, and 
4C) fell in the primary beam of the antenna, during 
the observation we tracked their geometric barycenter. 
As phase center s , we used the coordinates reported by 
Rodriguez et all (|1999l ) (Table [T]) in their VLA survey. 
Data were recorded in dual circular polarization with a 
velocity resolution (i.e. channel width) of ^0.1 kms -1 
and a total velocity coverage of about 55 kms -1 . To 
be able to measure and compare the relative position of 
the water maser emission across the epochs, we used the 
phase-reference technique using the close phase calibrator 
0333+321 (1.5 degree away). For this, we adopted a cycle 
time of 1 minute which left us about 20 seconds of inte- 
gration time on our calibrator and our target including 
slewing time of the antennas, for a total coherent integra- 
tion time on our target sources of about 2.5 h. 

Data were correlated at the VLBA correlator in 
Socorro (New Mexico, USA) we made three separate 
passes (one per source), using the coordinates of each 
three sources (IRAS 4 A,B,C) with an integration time 
of 2 seconds. The data reduction was performed using the 
Astronomical Image Processing System (AIPS) package. 
The determination of the phase calibration was done fol- 
lowing the standard method for spectral line VLBI data. 
The amplitude was calibrated using the template spec- 
tra method. The calibrator used to calculate the bandpass 
correction was 3C84, and the single-band delay corrections 
which was later applied to the spectral line sources was de- 
rived from the ICRF calibrator 0333+321. The fringe rates 
were estimate by selecting a bright and simple-structured 
channel and residual phase errors were removed later in a 
self-calibration process. The corrections were subsequently 
applied to the maser source to generate the channel maps. 
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Fig. 3. Spectra of the water maser emission in IRAS 4B as- 
sociated with the red-shifted northern spots and observed 
on June 11th and July 07th (see Fi gs. [1] and Et- The amb i- 
ent LSR velocity is +7.25 km s" 1 (jMardones et al.lll997l ). 
The peak velocity (km s^ 1 ) is indicated in each box. 
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Fig. 4. Spectra of the water maser emission in IRAS 4B as- 
sociated with the blue-shifted southern spots and observed 
on June 11th and July 07th (see Fi gs-fflandUl). The amb i- 
cnt LSR velocity is +7.25 km s" 1 (jMardones et al.lll997h . 
The peak velocity (km s -1 ) is indicated in each box. 



The final restoring beam of 1.2 x 1.2 mas and each channel 
was clean down to a final residual rms noise of about 15 
mJy. 



4. Results and discussion 

In the first epoch (Apr 1st), no emission was detected to- 
wards the four sources. In particular, IRAS 4C has never 
been detected in H2O, confirming the results from previ- 
ous interferometric studies (see Sect. 2). The other pro- 
tostars (IRAS4A1+4A2 and IRAS4B) presented water 
emission in two of the four epochs observed. This con- 
firms the high variability of the H qO spots associated 



with low-mass young stellar objects (IWilking et al. 1994; 



Claussen et all 11990 : Marvel et all 12008). with an active 
phase shorter than four weeks. Table [5] summarises the po- 
sition, radial velocity, and flux of the observed H 2 spots. 



4.1. NGC1333-IRAS4B 

In IRAS4B water emission was detected on June 11th 
and July 07th. Figure Q] reports the July 07th map: 
up to 13 H2O spots (labelled 'a' to 'm') showing a re- 
markable degree of collimation are seen. As a reference, 
the cross in Fig. [T] repor ts the coordinates of B follow- 
ing iReipurth et all (|2002j) who gave the position of the 
NGC 1333 protostars, observed at 3.6 cm, with a 50 mas 
uncertainty (dotted circle). The spot 'f is the closest one 
to the B coordinates, being inside the uncertainty area. 
Two boxes located at about (-70 mas, +100 mas) and (+70 
mas, -130 mas) mark the regions where the majority of 
the spots have been detected and that are shown in the 
zoom-in of Fig. [2l In particular, these spots ('a' to 'c': 
Upper panels; 'g' to 'm': Lower panels) arc the only ones 
detected in two epochs (Jun 11th and July 07th). The high 
degree of variability of the H2 O maser emission is clearly 
shown by the comparison between the integrated emission 
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maps of different epochs, as well as the H 2 line profiles 
reported in Figs. [3] and 2] . Note that the present angular 
resolution does not allow us to disentangle spots 'i' and 'j' 
as well as spots 'k' and '1'. However, these spots are well 
identified by different radial velocities, as shown by the 
spectra and by the channel maps of Fig. [5] 

The NW water emission (spots 'a' and T) is red- 
shifted, with velocities b etween +11.9 and +16 .8 km s _1 
(Vlsr = +7.25 km s" 1 ; iMardones et all (|l997l) . whereas 



the SE emission (spots 'g' and 'm') is blue-shifted (-0.1 
to +0.4 km s" 1 ). Thus, the positions and velocities of the 
water spots suggest that they are tracing a jet with a posi- 
tion angle of P.A. ~ -32° and driven by the 4B protostar. 

The pattern formed by the 'a', 'b', and 'c' spots has 
been identified in both Jun 11th and July 07th epochs. 
Consequently, we have been able to derive the apparent 
proper motion (see Fig. [6]), once having taken into account 
the parallax motion (including differential Galactic rota- 
tion) of such emission. Water masers are moving towards 
the NW direction with projected velocities in the range 
10-40 km s _1 (see black arrows). The SE spots show dif- 
ferent patterns in the two epochs: we can therefore only 
arbitrarily assume they are tracing the same components 
and then derive a typical proper motion: ~ 40-50 km 
s^ 1 . Table [3] reports the derived proper motions. These 
values are w e ll in a greement with the results found by 
" ( 2008L who found with their observations 



Marvel et al 



performed in 1998 a total separation velocity of 78 km 
s -1 . Note that for both the NW and SW spots, the mea- 
sured directions have a position angle between -30° and 
-45°, in agreement with the jet position angle, again in- 
dicating a very collimated flow from the B protostar. By 
using the measured proper motions and the radial velocity 
of the spots, we derive a high inclination to the plane of 
the sky between 10° and 35°, in agreement with the ob- 
served definite separation between the red and blue lobes 
of the jet. 

Figure [7] compares the position angle a nd the length 
of the H 2 jet as observed in 1998 by iMarvel et al 



(|2008h and in 2003 by us. Since the 1998 data are not 
in phase-reference mode we cannot positionally overlap 
the two maps and thus the plot of Fig. [7] is just indica- 
tive. However, it is possible to see that, after 5 years, 
the position angle, taking into account the corresponding 
uncertainties, remains nearly the same (-29° vs. -32°). 
Interestingly, a chain of water maser spots with the same 
inclina tion (-29°) was observed in 2006 by IPark & ChoJ 
(|2007l ). see their Fig. 3 in their VLA survey. 

The morphology and the kinematics of the H2O spots 
can be compared with the molecular outflow detected in 
several species (CO, CS, HCN, CH 3 OH, Si O H?CO) at 
larg er scales (HPBW s betw een 1" and 5") bv lChoil <|200lh 
and Jprgensen et al.l ( 2007 ) and thought to be driven by 
IRAS4B. The molecular maps indicate a clear bipolar 
structure, with a northern red-shifted lobe and a southern 
blue-shifted lobe, as shown by the H2O jet. The present 
maps suggest that, thanks to the H 2 spots, we are trac- 
ing the jet sweeping up the surrounding gas and creating 



the larger scale molecular outflow. However, the molec- 
ular outflow is located along the north-south direction, 
although the blue lobe slightly bends to S-SE. Therefore 
its position angle (~ 0°) is different from that of the H 2 
chain. As the inclination angle of 10-35 degrees excludes 
projection effects, either the molecular outflow is driven 
by another jet, not detected in H 2 0, and driven by an un- 
resolved companion of IRAS 4B, or the difference between 
the position angles is due to deflection or precession. In 
the first case, the jet could be bended by its interaction 
with the surrounding high-density medium. However, this 
effect seems to be ruled out for two main reasons: (i) the 
northern red-shifted and the southern blue-shifted lobes 
are rotated by t he same an g le, and (ii) according t o C 18 



and CS maps (jChoil l200lh iJorgensen et al.1 120071 ). there 



is no evidence of high-density clumps located along the 
main axis which could deflect the jet direction. On the 
other hand, if precession has occured the jet should have 
changed about 30° in around 110 years, which is the dy- 
namical scale of the molecula r outfl ow, using the arcsec- 
scale maps of ljprgensen et al] (|2007l) and assuming a typ- 
ical velocity of 100 km s _1 . In fact, the dynamical scale 
of the maser emission is of the order of a year. In other 
words, we should conclude that the IRAS 4B is processing 
wit h a definitely high r ate: ~ 18' yr _1 . 

Mar vel et al. (2008) noted that, in the case where 



the large scale outflow was created 110 years ago by 
a clock-wise precessing jet observed now in the water 
emission, the H 2 spots could be created by the im- 
pact of the jet against the walls of the outflow cavi- 
ties. In this scenario, the maser spots, whi ch require ex- 



trem e ly high-density (~ 10 -10 cm , e.g. lElitzur et al 



1989; Kaufman fc Neufeldl 1996) conditions, could orig 



nate in shocked layers compressed and accelerated at the 
interface with the surrounding medium. In particular, the 
maser amplification could be favoured at the border of 
the cavity, where the optical path along the line of sight 
is longer. This would also explain the high variability as- 
sociated with water masers. Actually, if we assume that 
maser emission is produced in the thin layer at the bor- 
der of the cavity, any small change of geometry associated 
with these regions, expected to be turbulent represent- 
ing the interface between mass loss and molecular cloud, 
would definitely change the optical path, thus affecting the 
maser intensity. This possible geometrical effect also would 
produce a change in the projected velocity, as confirmed 
by the spect r a rep orted in Figs. [3] and 2J Interestingly, 



Marvel et all (|2008l ) predicted that further VLBA obser- 



vations would reveal new masers south (north) of the NW 
(S E) maser groups. Th e comparison of our data with those 
of IMarvel et al. ( 2008 ) presented in Fig. [7] seems to sup- 
port this scenario, indicating for the H 2 map obtained in 
2003 a separation between the red- and blue-shifted em is- 
sion smaller than that observed by IMarvel et al. ( 20081 ) in 
1998. Future phase referencing VLBA observations would 
be instructive to check the validity of this hypothesis. 

Finally, it is tempting to give an estimate of the width 
of the H 2 jet by using the angle defined by the positions 
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Table 2. Offset position to the map center coordinates, 
radial velocities and peak fluxes of the observed H 2 spots 



100 - 



Spot 


epoch 


OJ2000 


f>J2000 




771 a 
F peak 




(2003) 


(mas) 


(mas) 


(km s^ 1 ) 


(mJy) 


IRAS 4B 


a 


Jun 11 


66.0 


95.0 


+15.6(0.1) 


21435(1850) 




Jul 07 


66.8 


96.0 


+15.9(0.1) 


981(6) 


b 


Jun 11 


65.8 


93.2 


+15.6(0.1) 


5305(45) 




Jul 07 


66.2 


93.8 


+15.1(0.2) 


2603(19) 


c 


Jun 11 


66.5 


97.0 


+15.8(0.1) 


2375(110) 




Jul 07 


68.4 


98.5 


+15.5(0.1) 


1709(12) 


d 


Jun 11 


- 


- 


- 


< 16 




Jul 07 


-33.0 


40.0 


+16.8(0.1) 


25013(1202) 


c 


Jun 11 


- 


- 


- 


< 8 




Jul 07 


-75.0 


120.0 


+17.0(0.1) 


850(59) 


f 


Jun 11 


- 


- 


- 


< 20 




Jul 07 


-28.0 


26.0 


+11.9(0.2) 


230(25) 


g 


Jun 11 


89.0 


135 


+0.2(0.2) 


31(2) 




Jul 07 








< 5 


h 


Jun 11 


89.1 


136.6 


-0.1(0.1) 


164(5) 




Jul 07 








< 4 


i 


Jun 11 


90.1 


137.5 


+0.2(0.1) 


40(2) 




Jul 07 








< 5 


j 


Jun 11 








< 6 




Jul 07 


90.5 


-138.0 


-0.2(0.1) 


40(2) 


k 


Jun 11 








< 3 




Jul 07 


90.7 


-138.5 


+0.1(0.2) 


49(4) 


1 


Jun 11 








< 6 




Jul 07 


91.55 


-138.5 


+0.4(0.2) 


13(2) 


in 


Jun 11 








< 3 




Jul 07 


91.55 


-139.2 


0.0(0.1) 


31(3) 


IRAS 4A2 


a 


May 14 


00.0 


00.0 


+8.3(0.1) 


330455(16010) 




Jun 11 


-0.2 


2.0 


+5.8(0.1) 


62545(4913) 


b 


May 14 








< 19 




Jun 11 


-3.1 


7.0 


+8.4(0.1) 


10512(990) 



a When H2O emission is detected, we report here the peak 
flux and the corresponding Gaussian fit uncertainty, 
otherwise we report the ler r.m.s. as upper limit. 

Table 3. Derived proper motions of the IRAS 4B maser 



spots 


Spot 


Aa 


AS 


v a 


Vt 




(mas) 


(mas) 


(km s" 1 ) 


(km s" 1 ) 


a 


-0.96 


+1.33 


15.1 W 


20.9 N 


b 


-0.59 


+0.53 


9.3 W 


0.8 N 


c 


-1.86 


+1.73 


28.9 W 


27.1 N 


h,k 


+1.50 


-3.02 


23.5 E 


47.0 S 


i,m 


+1.50 


-1.91 


23.5 E 


28.9 S 



of IRAS4B and spots 'e' and T (16°), and then measur- 
ing the width of this cone at the border of the IRAS 4B 
uncertainty area. In this way we obtain a width of 2 AU at 
a distance of ~ 12 AU from the driving protostar. Clearly, 
we cannot exclude that we are measuring the width of a 
cavity swept-up by the jet itself. However, as far as we 
know, this is the first attempt to estimate a Class jet 
width so close to the protostar. Another estimate has been 



o 
a 



-200 



-300 
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: IRAS4B 

- 
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* 
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1 
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■ 

t 

1 1 1 1 
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= Absolute star position 

■ 1998.Sep-Nov 

(Marvel et al., 2008] 
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300 200 100 
R.A. Offset (mas) 



-100 



Fig. 7. Comparison b etween the position a ngle of the H2O 
chain as observed by iMarvel et al. ( 2008h without phase 
referencing (black squares) and as reported in the present 
paper (2003 July; filled squares; red and blue in the elec- 
tronic edition only). After 5 years, the position angle re- 
mains the same (see text). The large circle is the position 
and corresponding uncertainty of the B protostar. 



inf erred for the HH212 je t by using the H2O spots imaged 
bv lClaussen et all (|l998h . In HH212 the width is 18 AU at 
50 AU from the driving source, indicating a decrease of the 
jet collimation moving away from the driving protostar. If 
we compare the IRAS 4B jet width with the range spanned 
by atomic microjets from more evolved Cl ass II sources 



surro unded only by thin disks (see Fig. 2 of ICabrit et al 
2007), there is no evidence of a higher jet collimation in 
Class sources embedded in ve ry dense cores. Thi s con- 
firms the conclusion reported bv lCabrit et al. ( 2007 ). who 
concluded that collimation is not due to external pres- 
sure, requiring that jets from young stellar objects are 
self-collimated by internal magnetic stresses. 

4.2. NGC1333-IRAS4A1+A2 

In IRAS 4A, maser emission was detected in the second 
and third epochs (May 14th, Jun 11th). Figured] (Left 
panel) summarizes the obtained maps: only two spots 
(named 'a' and 'b') were detected. The maps are centered 
on the coordinates of spot a as detected in May 14th: 012000 
= 03 h 29 m KK4T97221, <5 2 ooo = +31° 13' 32'.'286030. The 
velocities of the spots lie in the +5. 8, +6. 4 km s _1 range, 
thus being slightly blue-shi fted with respect to th e ambi- 
ent velocity (+7.25 km s _1 , iMardones et al. 1997). Figure 
|8] (Right panel) reports the highly variable H 2 spectra. 
Although the positions of the protostars are known with 
higher uncertainty (50 mas) than the H 2 spots (1 mas), 
the comparison between these positions shows that the 
maser spots are very close to the A2 protostar, being lo- 
cated at about 17 AU in the northern direction (offset: 
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+ 16 mas, -74 mas). On the other hand, Al is definitely 
more distant from the maser spots, at about 430 AU in 
the SE direction (+1405 mas, -1235 mas). This confirms 
the association o f mase r emission with A2, as also found 
by IPark fe Choil ( 20071 ). These authors interpreted their 
VLA maps as mainly a sign of a circumstellar disk with, 
in addition, a few spots related to a jet (see their Fig. 2). 
Unfortunately, despite the high resolution of the present 
data, since we have detected only two spots we cannot 
infer any information about the morphology and the kine- 
matics of the masering gas or discuss if the A2 water 
masers are tracing either a disk, a jet or a wide-angle cav- 
ity, as discussed for the IRAS 4B case. 

5. Summary 

We imaged the group of protostars in the NGC 1333 
IRAS 4 region in the H2O maser line at 22.2 GHz by using 
the VLBA in phase-reference mode. The maps allow us to 
detect a chain of water maser spots tracing the fast (> 
50 km s^ 1 ) jet driven by the IRAS4B protostar. We ob- 
served both the blue- and the red-shifted lobes very close 
to IRAS4B, < 35 AU. The H 2 chain is extremely colli- 
mated suggesting a jet width of ~ 2 AU. The comparison 
with the large scale molecular outflow suggests the jet is 
precessing with a rate of 18' yr . 
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Fig. 1. Water maser map from Jul 07th observations showing a chain of 13 H2O spots distributed over 73 AU, 
symmetrically locat ed with respect to th e driving protostar IRAS4B. The map is centered on the coordinates of 
IRAS4B fX-svmbol; lReipurth et ai1l2002l ): a 2 ooo = 03 h 29 m 12?003, 6 200 o = +31° 13' 08'.'14. The dotted circle indicates 
the uncertainty (50 mas) associated with the position of the protostar. The maser spots have been labelled with letters 
from 'a' to 'm'. For the sake of clarity, only spots 'd', 'e' and T have been identified. Two boxes identify the two 
regions where the rest of the spots have been detected, as shown by the zoom-in of Fig. [2 . First contour and step are 
~ 3 (0.84 Jy km s _1 ) and 6a, respectively. 
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Fig. 2. Zoom-in of the two regions identified by black boxes in the Jul 07th map of Fig. Q] . The May 14th (no 
detection), June 11th, and July 07th epochs are reported. The la level is 0.15 Jy km s _1 for the June 11th map, and 
0.28 Jy km s _1 for the July 07th map. Upper panels: Northern region associated with the red-shifted (see Fig. [3]) spots 
labelled 'a', 'b', and 'c'. First levels and steps are 3 and 12a, respectively. Lower panels: Southern region associated 
with the blue-shifted (see Fig. 2]) spots labelled 'g', 'h', 'i', 'j', 'k', T, and 'm'. First levels and steps are 3 and 6a, 
respectively. Note that spots 'i' and 'j', as well as spots 'k' and '1', are not well separated in the present integrated 
H2O emission map, but they are associated with different velocities (see channel maps of Fig. [5]). 



Desmurs et al.: The H 2 jet from the NGC 1333-IRAS4B protostar 



9 



RAS4B South 



-136 



z~i 1 1 1 1 - 



-138 - 



-140 
■136 



-138 



-140 



-0.51 



I i 1 1 1 r~ 



-0.40 




=H 1 1 1 \-^ = 



+0.02 



I i 1 1 1 r~ 



-0.29 




■H h 



H h^=H h 



+ 0.13 




I I 1 1 1 1 - 



-0.19 




H h^=i 1 



+ 0.23 




I I 1 1 1 1 — 



-0.08 




H h^=H h 



+ 0.44 




J I I I I I I I I I I I L 




H \-= 



+0.65 



91.5 90 91.5 90 91.5 90 91.5 90 91.5 

R.A. Offset (mas) 



90 



Fig. 5. Channel maps of the H2O emission towards the IRAS4B southern region, where the blue-shifted spots have 
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Fig. 6. Zoom-in of the northern (Left: spots 'a', 'b', and 'c') and southern (Right: spots T to '1') IRAS4B regions, 
where H2O emission has been detected. The contour maps show the comparison between the June 11th (dotted line) 
and July 07th (solid line) epochs. For the sake of clarity only selected contour levels have been drawn. Note that spots 
'i' and 'j' as well as spots 'k' and T are not fully spatially resolved, but are associated with different velocities (see 
spectra in Fig. [4] and channel maps of Fig. [5]) . Solid black arrows indicate the apparent motions (parallax-subtracted) : 
the water spots are moving away from the protostar at projected velocities of between 10 and 52 km s^ 1 . 
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Fig. 8. Left panel: Water maser maps measured for IRAS4A (dotted contour: May 14th; solid contour: Jun 11th). 
First contour and step are ~ 3 (1.5 Jy km s^ 1 ) and 13a, respectively. The A2 protostar (02000 = 03 h 29 m 10M21, 
^2000 = +31° 13' 32'.' 21) is located at +16,-74 mas, whereas the Al protostar is at +1405, -1235 mas, according to 
Reipurth et"aH 1)20021 ). Right panels: Spectra of the H2O maser spots: the pea k velocity (km s 1 ) is reported in each 



box. The vertical dashed line is for the ambient LSR velocity (+7.25 km s ; iMardones et al.lll997h 



